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Abstract 

It has been suggested by Bartlett and Hall that our universe may have the 
critical density in baryons by virtue of specific interactions with a 'shadow' 
world. We show that this possibility is severely constrained by primordial 
nucleosynthesis, stellar evolution and the thermalization of the cosmic mi- 
crowave background. In particular, recent observations of small angular-scale 
anisotropy in the cosmic microwave background conclusively rule out all such 
baryon-dominated cosmologies. 
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I. INTRODUCTION 



The 'fine tuning' problems of the Friedmann-Robertson- Walker cosmology favour the possibility 
that the universe is spatially flat with density parameter O = 1, as would result naturally from 
an inflationary phase in the early universe However, the observed abundance of deuterium, 
which could only have been formed during big bang nucleosynthesis (BBN) [j^], implies an upper 
bound of f^N < 0.033/i~ 2 on the nucleonic contribution to the density parameter ||. Since the 
Hubble parameter h (= -ffo/100 kmsec _1 Mpc _1 ) is in the range ~ 0.5 — 1, ordinary matter cannot 
close the universe and this discrepancy has motivated the much discussed idea that non-baryonic 
dark matter makes up the difference.^] This would require the existence of new stable massive 
particles in extensions of physics beyond the Standard Model, e.g. massive neutrinos or the lightest 
supersymmetric particle ||. It is clearly important to establish whether there are any loopholes in 
this argument, especially since nucleons can in principle make up all the dynamically detected dark 
matter which amounts to only Q ~ 0.1 |1],[<J|. In particular, recent observations of microlensing of 
stars in the Large Magellanic Cloud indicate that ~ 0.5 Mq objects (probably old white dwarfs) 
can already account for upto half of the dark matter in the halo of our Galaxy Q . 

There have been several proposals for relaxing the BBN bound on the nucleon density by either 
modifying general relativity ||, or complicating the standard picture of nucleosynthesis, e.g. by 
invoking non-zero neutrino chemical potentials or inhomogeneities in the nucleon distribution 
resulting from a first-order quark-hadron phase transition [11]. The latter possibility is physically 
the best motivated but detailed studies show that even with new free parameters it is difficult to 
permit the critical density in nucleons [12]. There have also been suggestions which invoke physics 
beyond the Standard Model fl3| . One requires a tau neutrino with a mass of 20-30 MeV (but with 
its relic abundance suppressed by a factor of ~ 10 below the nominal 'freeze-out' value) decaying 
'invisibly' into electron neutrinos with a lifetime of 200-1000 sec [14|. Although this would allow 
a = 1 universe, the model is rather fine-tuned and quite unlikely from the particle physics 
viewpoint; in particular the v T must have its tree-level Standard Model decays {y T — > e + e~v e ) 
totally suppressed as these would disrupt nucleosynthesis [15,13]; moreover the required mass 
region is now almost ruled out experimentally |tl6| . A better motivated suggestion [jl^| is that the 
hadronic and radiation cascades triggered by the decay after BBN of massive unstable particles such 
as gravitinos can process the elemental yields so as to synthesize adequate amounts of deuterium 
even for = 1. However, there is a problem with the abundance of 6 Li, which, in this model, is 
nearly as abundant as 7 Li, in apparent conflict with observations |l8|]. (The same comment applies 
to the model invoking Hawking radiation from primordial black holes [19].) 

In this paper we consider another exotic idea due to Bartlett and Hall [^(J viz. that our universe 
(the 'visible' sector) may have become coupled to a 'hidden' sector through a phase transition after 
nucleosynthesis. The mixing of photons in the visible sector with their counterpart 'paraphotons' 
in the hidden sector can then result in a sudden decrease of the comoving photon density. Thus 
the universe may well have the critical density in nucleons today, while having a nucleon-to-photon 
ratio during BBN consistent with observed elemental abundances. To see this we note that the 



There may exist new forms of baryonic matter which do not participate in nucleosynthesis, e.g. strange 
quark matter 'nuggets' Q. If these survived baryon evaporation, they could make a substantial contribution 
to £1 today & therefore we distinguish between baryons and nucleons in the context of BBN. 
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nucleon mass density in ratio to the critical density is 

__ m N n N _ 87rm N f n N \ 2((3)T$ 
UN ~ Pc -3H$M*{n 7 ) ' W 

where the present temperature of the blackbody cosmic microwave background (CMB) is To = 
2.728 ± 0.004 K (95% c.l.) [0. Therefore we require 

~ 2.72 x W- 8 h 2 e- 3 , (2) 

Q N =1 

for a critical density in nucleons, where G = To/2.73K. In the standard cosmology, the comoving 
entropy does not change after e + e~ annihilation so the nucleon-to-photon ratio during nucleosyn- 
thesis equals its value today, i.e. t/bbn = ?7o- As mentioned earlier, the observational lower limit 
to the abundance of deuterium (together with observational upper limits to the abundances of 
helium and lithium) requires ?7bbn < 8.9 x 10~ 10 §. The interesting possibility suggested in ref. 
p0|| is that the comoving photon number decreases following nucleosynthesis, so that rjo is raised 
to ~ (10 — 100) t?bbn and can be as high as rj c . 

We examine the effects of such physics on various cosmological and astrophysical processes such 
as BBN itself, the spectrum of the CMB and stellar evolution. It turns out that these arguments 
severely constrain but cannot quite rule out the scenario. However, all models in which the critical 
density is comprised entirely by baryons (nucleonic or otherwise) are shown to be definitively 
excluded by recent measurements of anisotropy in the CMB together with measurements of the 
power spectrum of large-scale structure (LSS) in the universe. Thus, there appears to be no 
alternative to non-baryonic dark matter if indeed f2 > 0.1. 
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II. THE PHOTON COOLING COSMOLOGY 

First we briefly review the scenario of Bartlett and Hall fl20| . They consider a gauge theory 
with the group G®G' at energies far above the electroweak scale (v ~ 246 GeV). The group G 
contains the Standard SU (3) 03 SU (2) (g) U(l) Model while G' corresponds to a hidden sector with 
no tree-level renormalizable couplings to known particles although similar in its symmetry breaking 
pattern. In particular G' breaks to a group which includes a hidden (unbroken) electromagnetism 
U(l)' at a scale v' <C v and the 0(100) light degrees of freedom in G' (generically termed X') 
acquire a mass at this scale. The cosmology is thus quite novel. The lightest 'connector' particle 
(carrying both G and G' quantum numbers) has a mass Mq > i; so that scattering processes 
involving the connector particles maintain the two sectors in thermal equilibrium at T > Mq. 
At lower temperatures these reactions freeze-out and the two sectors decouple. Subsequent phase 
transitions and annihilations of heavy species affect the sectors differently resulting in T' < T. 
(For example if the temperature difference is dominantly due to heavy particle annihilations in 
the visible sector then T/T' ~ 2 (11/4) x / 3 ~ 3 at v' < T -C m e @.) Thus, during BBN the 
universe contains two plasmas which are not in thermal contact with each other. At some point 
after nucleosynthesis a 'recoupling' reaction reestablishes thermal contact, cooling the photons 
in the visible sector and hence reducing the comoving photon number density. Such a reaction 
requires a renormalizable interaction between photons and the X' particles in the hidden sector. 
Indeed the particles X' can have a very small electric charge, which may arise naturally in a unified 
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non-Abelian gauge theory [22], so that the Compton scattering reaction 7 J' ^ ^ X 1 shown in 
Figure 1 can recouple the two sectors. 

The rate of this process may be estimated as follows. Wave function mixing mediated through a 
one-loop diagram involving the connector particles generates a term {e/2)F txv F' in the Lagrangian, 
where e ~ (ee' /I6ir 2 )vv' /Mq for Mq 3> v,v' and values of e ^ 10~ 3 may be considered natural 



[20]. After diagonalizing the kinetic terms, the interactions of the photons can be written in several 
alternative bases. Let and J ,fl be the currents, including the gauge coupling constants e and 
e', for U(l) and U(l)' respectively. When only ordinary matter is relevant it is most convenient to 
write the corresponding Lagrangian to order e in the form 

Ant = a% + (a'" + ea^) (3) 

where is defined to be the field coupled to which creates "our" photon 7. (Note that the 
quanta coupled to and J 1 * 1 are not orthogonal. Thus the field a'^ does not create the shadow 
photon 7' of the hidden sector; this is done instead by a field A /fl .) From eq.(||) we see that the 
hidden particles X' interact with the photon through a very small electromagnetic charge e e' . For 



fermionic X' species the rate of the recoupling process is then [20] 



n x , (av) = C n x > (4) 

where a' = e' 2 /Air and C is a number of 0(1) which comes from thermally averaging over the 
energy transfer. As we shall see, constraints from the CMB require the recoupling reaction to 
occur during the radiation-dominated era. The energy density of the universe is then dominated 
by the plasma in the visible sector. Defining the recoupling temperature Tr to be the temperature 
of our plasma when the recoupling rate equals the expansion rate gives 



where a = e 2 a' X]* g[ if and g[ is the number of spin states (assumed fermionic) of the hidden plasma 
having interaction qi e' with 7'. To calculate the resultant cooling of the visible sector it is assumed 
that the energy transfer occurs quickly compared with the expansion rate; the final temperature 
Tp is then obtained by demanding energy conservation at recoupling, i.e. gT^ + g'T^ = (g + g')T^. 
The energy density of the shadow sector should be negligible during BBN (in order not to increase 
the expansion rate excessively) but dominates after recoupling (because g' S> g = 2) which gives 
n 7 (Tp)/n 7 (7R,) ~ (g/g') 3 ^ 4 - Therefore a critical density nucleonic universe (see eq. (|2j) ) requires 
g' ~ 40 — 1000 and, correspondingly, T/T' > 2 — 5. Subsequently the X' particles acquire mass 
through a phase transition with a critical temperature bounded as 

< 2r < T R /2 . (6) 
It would be natural to expect that mx 1 is of 0(T' C ). 

A. Constraint from BBN 

Nucleosynthesis obviously provides an upper bound on the recoupling temperature since the 
abundances of the elements would be altered if it is too high. To study this, we have modified the 
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standard computer code [23| to include an abrupt drop in temperature at Tr when the recoupling 
reaction occurs. Before this, the universe contains two plasmas which are not in thermal contact, 
the usual one characterized by the photon temperature T = T 7 (consisting of 7, 3 v's, e^, as 
well as nucleons), and the other one containing the shadow sector particles like 7' and X' with 
temperature T". An important parameter is the ratio of the shadow sector temperature to the 
neutrino temperature before recoupling, 

T' 

r^Tfr, (7) 

which stays constant during e + e~ annihilation when the photon temperature T increases (but not 
T' or T u ). It enables the energy density p' of the shadow sector before recoupling to be tracked 
using the neutrino energy density, p' = p u g'r 4 Jg v , and is thus a more convenient parameter than 
T'/T as used in ref. f20|| . Now at the recoupling temperature Tr there is a sudden drop (which 
would be smooth in a more detailed model) to the final temperature Tp. Energy conservation gives 

fi(<7 7 + 9u + 5c)T R + //(T R ) = jg(s 7 + g v + g e )Tp + ^g% , (8) 

where g-y = 2, g u = 21/4 and g e = 7/2. Since Tr = Tp for neutrinos, the corresponding terms above 
cancel and we obtain (neglecting the electron energy density for the temperatures of relevance), 



T F = AT R , A 



2 + (A)4/3 gV 4 nV4 

2 + g' ■ W 



After recoupling the energy density of the shadow sector is just p' = ^g'T^ since now T = T' . 

We have implemented a routine in the BBN code dividing the thermal history into the three 
steps (before, during and after recoupling) and compute the expansion rate accordingly, taking 
into account the contribution of shadow particles. We assume 3 light neutrino species, use the 
neutron lifetime r ra = 887 db 2 s [jT6[| , and incorporate small corrections to the helium abundance as 
reviewed in ref. . Since the main effect of photon cooling is to increase the helium-4 abundance 
and decrease the deuterium abundance, we bound their values from recent observations of helium 



in metal-poor extragalactic HII regions [24] and of deuterium in the interstellar medium [25|: 



y P ( 4 He) < 0.25, (D/H) p > 1.1 x 1(T 5 . (10) 

We do not consider the bound on 7 Li because it does not provide a useful constraint. The adopted 
bounds are very conservative, based upon consideration of a variety of data which are critically 
discussed elsewhere fi"3|| . As shown in Figure 2, we can then derive the maximum value of Tr 
allowed for each given set of parameters r and g' . The allowed region is: 

T R < 74 keV , r < 0.3 and 40 < g' < 2340 , (11) 

only part of which is shown in the figure for convenience. (The highest allowed value of g' corre- 
sponds to the smallest initial value of rj, viz. t/bbn > 1-7 X 10~ 10 [H].) For comparison, the values 



suggested in ref. §0| were T R < 30 keV, r < 0.7 and g' ~ 40 - 1000. 
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B. Constraint from the CMB spectrum 



A lower bound on the recoupling temperature follows from the observation by COBE that 



the CMB spectrum is accurately blackbody [21|. The recoupling reaction removes photons from 
the thermal bath thus inducing a chemical potential for the photon distribution. If the reaction 
occurs early enough, however, photon creation processes can restore the Planck spectrum |26| in 
agreement with the observational upper limit [21| on such a chemical potential: /j<9x 1CP 5 (95% 



c.l.). To thus constrain Tr we use a simple argument concerning the relevant time scales. 

At temperatures below the e + e~ threshold, the processes which determine the spectral shape 
are Compton scattering (C) which preserves the photon number, and the radiative processes 
bremsstrahlung (B) and double Compton scattering (DC). The time scales of these processes are ^ 



t c = 9.25 x 10 w sec 



1 



t D c = 2.53 x 10" 2 sec - 



rj \keV 



1 



t B = 2.18 x l(T 8 sec -=■ 



rj \keV 



1 



rj 2 \keV 



-5/2 



In 



-2 



1 / T y/ 4 
^ UeV 



(12) 



For a On = 1 universe, the value of rj after recoupling (i.e. today) is given by eq. (|2j) whereas before 
recoupling it is 



rj{T > Tr) = rjo 



rj Q A 6 



(13) 



where A is defined in eq.(^). (Note that baryons and nucleons are equivalent with regard to the 
CMB since only their electromagnetic interactions are relevant; therefore we henceforth refer to Or 
rather than On-) The time scales above are to be compared with the age of the universe obtained 
by integrating the Friedmann equation, 



2.42 x 10 6 sec g\ 
2.42 x 10 6 sec 



1/2 (_t_\ 

VkeV J 



+ 



Jd¥ (kev) 



for T > T R , 
for T <T F . 



(14) 



where we have used 



9i = 9 (T>T R )=2 + f (^) 4/3 + 5 V(A) 4/3 
gF ^ g (T<T F )=2 + fA^(±y /3 + g ' . 



(15) 



Obviously the abrupt transition at Tr would be smoothed out in a more detailed picture of the 
phase transition. 

In Figure || we show the various time scales for a typical set of parameters, with a sudden jump 
in rj as the temperature drops from Tr to T F as determined by eq.([|). Although bremsstrahlung 
usually dominates at low redshifts for a high density baryonic universe, the most efficient photon 



2 We have inserted the correction factor of An |2(| into the rate for bremsstrahlung [^7j . 
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creation process is always found to be double Compton scattering for the values of Tr, r and 
g' under consideration. The limiting value of the recoupling temperature is then conservatively 
obtained by demanding that the time scale for this be shorter than the expansion time scale. We 
have indicated in the figure the error bands for the time scales in eq.(|l2|) corresponding to the 
allowed region of rj after nucleosynthesis, viz. r\ G (1.7 — 8.9) x 1CP 10 Q. To obtain the minimum 
recoupling temperature we took the highest possible value for r/o from eq.@ and searched for the 
temperature at which t > £dc is satisfied as a function of the parameters r and </; this gives 

T R >0.73keV (^f) 1/3 • (16) 

Combining the constraints from BBN and the CMB spectrum it follows that the recoupling tem- 
perature must be in the range 

3keV < T R < 75keV (17) 
for the allowed range of r and g 1 as given in eq. (|TT[) . 



C. Constraints from stellar evolution 

Further constraints on the model follow from the energy loss argument for stars |[20|| . For 
example if the X' particles are light enough to be produced in red giants through plasmon decays, 
a limit on a (eq.(|5|)) follows from requiring that the energy loss through X' emission does not affect 



the stellar luminosity or age to a noticeable degree [28]. The resultant constraint is very stringent, 



viz. a 1CP 28 [ 20 1 , and in conflict with the required value from eq.(|5[) unless the recoupling 
temperature is very low. This, however, is not allowed by the CMB spectrum constraint (eq.(|l6|) ) 
so the only escape route is to argue that the X' particles are too heavy to be produced in stars |2(| . 
In the cores of red giants the plasmon mass is uj p \ ~ 20 keV, so this requires mx' > ^pi/2 ^ 10 keV. 
There are additional contributions to X' production from Compton and bremsstrahlung processes 
and to suppress these one would require my > 5T where T ~ 10 keV is the core temperature p(|.p| 
However, if my is indeed of 0(Tr) then this implies Tr > 50 keV (see eq.(^)), which exceeds the 
upper limit from BBN over most of the parameter space (see Figure Q). Thus it would be necessary 
to invoke some fine-tuning and argue that mx> 3> Tr. 

An even more restrictive constraint on a comes from supernova physics. After collapse, the 
supernova core forms a neutron star which radiates its binding energy primarily in the form of 
neutrinos which diffuse out of the dense central region over a time scale of seconds. Observations 
of (anti) neutrinos from SN1987A [29] confirmed that these indeed carried away ~ 10 53 erg on the 



expected time scale. A competing mechanism for cooling the core would be the production of X' 
particles through plasmon decays. Then a conservative requirement is that the X' particles do not 
carry away more energy than the neutrinos and thus curtail the duration and energetics of the 
neutrino burst. Assuming the core radius to be ~ 10 km and a cooling time of ~ 10 s one obtains 
for the energy loss E' through X' [20 



3 Note that processes of the form e 7 — > e 7' or e 4 He — > e 4 He7' are not allowed by the Lagrangian of 
eq.(|^); if these existed, the associated energy loss could not be evaded by increasing mx> ■ 
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E' « 10 71 erg affil% < 10 53 erg , (18) 

where T 10 = T/(10 MeV), p u = p/(W u g/cm 3 ), and the X' particles are assumed to be lighter 
than the plasmon mass in the supernova core. Typical values for the density and temperature are 
P14 ~ 6 and T\q ~ 5 |3(J. Solving this equation for a and inserting into eq.(||) we then find: 

T R < 3.6 x 10 5 keV p^ /3 T w 3 Ca' r 2 . (19) 

In Figure || we show that this constraint together with those from BBN and the CMB spectrum 
effectively rules out the model (taking C = 1, and a' = 0.1). The only allowed region is a tiny 
corner at the highest value of r and the lowest value of g' , fixing the recoupling temperature to be 
7-8 keV. This region is further limited by decreasing a' or C. Of course the physical parameters 
of supernova collapse are somewhat uncertain so this argument is not definitive. One can only 



conclude that X' particles, if they exist, may play an essential role in the supernova process [20]. 



D. Constraints from LSS and CMB anisotropy 



Finally we consider the formation of large scale structure and the concomitant generation of 
small scale CMB anisotropy in a critical density baryonic universe.^ In the standard picture, such 
structure grows by gravitational instability from a (nearly) scale-invariant spectrum of primordial 
adiabatic, gaussian fluctuations (with power spectrum P(k) oc k n , n ~ 1), presumably generated 
during an inflationary phase in the early universe On scales smaller than the horizon size 

at matter-radiation equality the primordial spectrum is modified and its slope becomes negative. 
The COBE detection of anisotropy in the CMB on angular scales exceeding the horizon size at the 
recombination epoch (3^] has provided strong support for this basic picture and moreover fixed the 
normalization of the primordial spectrum. Subsequent observations of anisotropy on small angular 
scales have also provided evidence for a 'Doppler peak' in the angular power spectrum such as 
would be expected from acoustic oscillations of the coupled plasma and photon fluids during the 
last scattering of the CMB J3 1[J - Besides the observed power spectrum of galaxy clustering, we thus 
have a new quantitative test for models of structure formation. 

A baryonic universe was of course the first to be investigated in the context of structure for- 
mation [^]. The salient feature of such a universe is that fluctuations can begin to grow only after 
(re)combination (rather than as soon as the universe becomes matter-dominated) since the baryons 
are tightly coupled to the photons before this epoch. Moreover the coupling is not perfect (and be- 
comes weaker as the universe turns neutral) so adiabatic fluctuations suffer 'Silk damping' on scales 
smaller than the damping mass of M D » 1.3 x lO 12 (ft B /i 2 )" 3/2 M . It was recognized early on that 
primordial gaussian fluctuations with spectral index n < 2 would thus result in excessively large 
fluctuations in the mass distribution on small scales []34j]. Moreover the upper bounds on CMB 
anisotropy available a decade ago already implied that there would not be sufficient time after 



4 Note that in the Bartlett-Hall model [^[J the X' particles of the shadow sector can be a candidate for 
the dark matter. Since nx> ~ ^7 at recoupling, there will be a large contribution to the density parameter 
(fix 1 ~ AQh~ 2 {mx' /keV)) unless the X' particles subsequently self-annihilate efficiently after becoming 
non-relativistic. This is however quite possible and, in fact, very likely in analogy to the 'visible' sector. 
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(re)combination for the implied small primordial fluctuations to grow to the large-scale structure 



observed today in a baryonic universe [ 35fl . 

In the past particle physicists have tended to ignore such constraints since there did not exist 
a 'Standard Model' for the formation of structure. The increasing weight of observational evidence 
for the picture outlined above, however, now makes such a viewpoint untenable. In particular the 
data (as opposed to upper limits) on CMB anisotropy now allow a precise quantitative test of a 
baryon-dominated universe. For a Ob = 1 universe, the slope of the primordial spectrum n and 
the Hubble parameter h are the only free parameters. Thus we will first determine the values of n 
and h which would allow a fit to the data on small-scale CMB anisotropy. These are available as 
measurements of the multipole moments C\ = J2 
the celestial sphere in spherical harmonics: 



m af m , expanding the CMB temperature on 



AT = ^2ai m Yi m (6, 



(20) 



Lm 



Then we can check if the power spectrum of galaxy clustering can be reproduced for the same 
choice of parameters. We emphasize that this test is thus as model-independent as possible. 

To compute the angular power spectrum for each model specified by a pair (n,/i), we use a 
fast Boltzmann code and adopt the COBE- normalization at large angular scales. Then we 
convolve each power spectrum with the experimental 'window functions' before comparison with 
the observational data by generating a x 2 surface over the (n, h) space under investigation; this 
method is outlined in detail in ref. [37]. The x 2 value is computed according to 



X 2 (n,h)= £ 
JV=1 



ST, data (N) 



-, 2 



a 



data 



(AO 



(21) 



where the sum is over all the CMB observations, 5T z model (iV, n, h) and 5T^ ta (N) are, respectively, 
the flat band powers (i.e. convolved with the experiment-specific 'window functions') for the models 
and observations, and <7 data (iV) are the associated observational errors. The CMB data we consider 
and the corresponding errors are tabulated in j3?| and shown in Figure ||. (We consider only the 
central value of the Saskatoon data which have a ±14% calibration uncertainty.) The band powers 
are given by 



I(Wi{N)) \ 



4-7T 



1/2 



(22) 



where the deconvolving factors I(Wi(N)) are the logarithmic integrals of the experimental window 
functions Wi{N) defined as 



J(Wi) = 2ir £ 



1=2 



(21 + l)Wi 
4ttI(1 + 1) ' 



(23) 



and l max can be safely taken to be 1200 for all existing data (for larger I the acceptance is negligible). 
Finally, the angular scale investigated by experiment iV is given according to 



les(N) 



IjlWijN)) 
I(Wi(N)) 



(24) 
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To obtain a x 2 value for each point of the (n,h) space we use the data and eq.(|^) for each 
observation N and sum over all existing observations. The values of the multipoles C[(n,h) are 
obtained by running the fast Boltzmann code mentioned above|] 

In Figure |6] we show the region of (n, h) space which is consistent with the data. We investigate 
the parameter range 0.25 < h < 1 and 0.5 < n < 1.5, the former covering all direct measurements 



[16] and the latter based on general theoretical considerations of inflationary models [38|. We show 



both the 68% C.L. and 95% C.L. contours; the latter delineate the range 

0.59 < n < 0.78 , 0.41 < h < 0.62 . (25) 

The negative correlation of n and h arises because the amplitude of the acoustic peak plays the 
dominant role in the fitting procedure and for J1b = 1 this peak amplitude decreases as h and 
n decrease. Thus a decrease in n accompanied by a correlated increase in h (or vice versa) can 
preserve the amplitude and thus the \ 2 value for the fit. In general rather low values of both n 
and h are required to fit the acoustic peak in a critical density baryonic universe. 

Given these favoured values of n and h we can now calculate the linear matter power spectrum. 
In Figure we show the spectra for the two extreme models in the (n, h) plane, viz. n = 0.5, 
h = 0.62 and n = 0.78, h = 0.25, along with the data points inferred from the APM galaxy 



survey [39|. To compute these spectra we have used the COSMICS code [4C] and normalized to 



COBE using the prescription given in ref. [41|. The strong oscillations seen at short scales arise 



from matter fluctuations on wavelengths smaller than the Hubble radius at (re) combination and 



reflect the value of the phase of the perturbation at recombination [34]. It is seen that even for 
extreme choices of parameters, the power spectra of J7b = 1 models disagree strongly with the 
observational data. Even though the clustering of APM galaxies may be enhanced over that of 
the matter distribution, the discrepancy cannot be resolved for any sensible value of the 'bias' 
parameter. The situation is further aggravated if the value of h is decreased in Figure ^(a) or the 
value of n is decreased in Figure |7|(b), staying within the allowed region of eq.(^), since in both 
cases the matter power spectrum is further damped at short scales. Thus we can firmly establish 
that it is not possible to simultaneously satisfy the LSS and CMB data in an universe with the 
critical density in baryons. 



III. DISCUSSION 

Although big bang nucleosynthesis is usually cited as providing the main argument for non- 
baryonic dark matter, it is possible to entertain deviations from the standard picture in order 
to accommodate a baryon-dominated critical density universe. One such suggestion invoking a 



reduction in the comoving entropy after nucleosynthesis [ 20 1 is particularly interesting because it 
challenges the standard assumption that the thermodynamic history of the universe is well known 
at least up until the BBN epoch. We have demonstrated that this model is very tightly constrained 



We ignore a possible contribution from gravitational waves (which would lower the COBE normalization) 
and adopt for the helium fraction the central value lp( 4 He) = 0.24 ]24[1 . Note also that the error bars 
cr data (./V) are mostly asymmetric; to accommodate this in the % 2 calculation, we establish for each data 
point whether the model lies above or below and then take the appropriate error bar in eq. (^2|) . 
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(if not ruled out altogether) from consideration of other processes such as thermalization of the 
CMB spectrum and the effects on stellar evolution of the new physics invoked by the model. Finally 
we have shown that recent observations of small-scale CMB anisotropy and large-scale structure 
decisively rule out this model. This conclusion also holds for other models which attempt to evade 
the constraint from nucleosynthesis in order to allow a baryonic Einstein-De Sitter universe. 
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FIGURES 




FIG. 1. The Bartlett-Hall mechanism for recoupling of the two sectors via photon mixing. 




FIG. 2. Constraints on the recoupling temperature as a function of the model parameters r 
and g' from BBN and the CMB spectrum. The allowed region is between the two surfaces. 
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FIG. 3. Time scales for Compton scattering, double Compton scattering and bremsstrahlung 
compared with the expansion time scale in the Bartlett-Hall model. The steps correspond to the 
recoupling process at which the temperature drops from Tr to Tp and the comoving photon density- 
decreases by a factor of ~ 10 — 100. 
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FIG. 4. Constraints on the recoupling temperature as a function of the model parameters r 
and g' from the CMB spectrum and supernova physics (for p\± = 6, T\q = 5, C = 1 and a' = 0.1). 
The region between the surfaces is excluded. The only allowed region obtains for the highest values 
of r and the lowest for </, fixing the recoupling temperature to be ~ 7 keV. 
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FIG. 5. An illustrative fit to the data on the CMB angular power spectrum of a S7b = 1 universe 
with primordial spectrum slope n = 0.66 and Hubble parameter h = 0.3). 
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FIG. 6. x 2 contours (68% and 95% C.L.) in the n — Hq plane for the fit to observations of CMB 
anisotropy of a S7b = 1 universe. The negative correlation between the primordial spectral index 
and the Hubble parameter is discussed in the text. 
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FIG. 7. Predicted power spectra of density fluctuations normalized to COBE for two Ob = 1 
models, (a) n = 0.5, h = 0.62 and (b) n = 0.78, h = 0.25, compared with APM data. 
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